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Abstract 


The y-ray emitting compact symmetric objects (CSOs) PKS 1718-649, NGC 3894, and TXS 0128+554 are lobe- 
dominated in the radio emission. In order to investigate their y-ray radiation properties, we analyze the ~14 yr 
Fermi/LAT observation data of the three CSOs. They all show the low luminosity (10*1-10? erg s~') and no 
significant variability in the y-ray band. Their y-ray average spectra can be well fitted by a power-law function. 
These properties of ^-rays are clearly different from the y-ray emitting CSOs CTD 135 and PKS 1413+135, for 
which the 7-rays are produced by a restarted aligned jet. In the L, — I’, plane, the three CSOs are also located at the 
region occupied by radio galaxies (RGs) while CTD 135 and PKS 1413+135 display a similar feature to blazars. 
Together with a similar radio emission property to y-ray emitting RGs Cen A and Fornax A, we speculate that the 
7y-rays of the three CSOs stem from their extended mini-lobes. The broadband spectral energy distributions of the 
three CSOs can be well explained by the two-zone leptonic model, where their y-rays are produced by the inverse 
Compton process of the relativistic electrons in extended regions. By extrapolating the observed Fermi/LAT 
spectra to the very high energy band, we find that TXS 01284-554 among the three CSOs may be detected by the 


Cherenkov Telescope Array in the future. 
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1. Introduction 


Compact symmetric objects (CSOs), a sub-class of active 
galactic nuclei (AGNs), are distinguished by a symmetric 
and compact radio structure with the projection size 
<1 kiloparsec (kpc). Their symmetric radio structure is thought 
to be due to a misaligned jet with no or weak relativistic effect 
(Phillips & Mutel 1980; Wilkinson et al. 1994; Readhead et al. 
1996). Thus CSOs look like the mini-version of off-axially 
observed radio galaxies (RGs). However, some CSOs show 
obvious flux variation at the radio band, even significant 
variability and high luminosity in y-ray band, indicating the 
strong Doppler boosting effect of a relativistic jet (Gan et al. 
2021, 2022). Recently, Kiehlmann et al. (2023) reported that in 
addition to the compact symmetric structure, no significant flux 
variation and superluminal motion should be taken as a criterion 
for the CSO selection, that is to say no relativistic effect. 

It is well known that blazars are the main one of the 
extragalactic y-ray emitters since a relativistic jet dominates their 
multi-band radiations, and their y-rays are from an aligned 
pc/sub-pc jet. Generally, it is believed that the 5-rays of RGs are 
also generated from their pc/sub-pc jets (e.g., Abramowski et al. 
2012; Aliu et al. 2012; Grandi et al. 2012; Casadio et al. 2015; 
Tanaka et al. 2015; Xue et al. 2017). Nevertheless, the detections 
of y-rays from the large-scale lobes of RGs Cen A and Fornax A 
by the Fermi/LAT indicate the existence of high-energy particles 


in these extended regions (Abdo et al. 2010; Ackermann et al. 
2016). Until now, five CSOs have been detected in the ^-ray 
band, including PKS 1718-649 (Migliori et al. 2016), NGC 3894 
(Principe et al. 2020), TXS 01284-554 (Lister et al. 2020), CTD 
135 (Gan et al. 2021), and PKS 14134-135 (Principe et al. 2021; 
Gan et al. 2022). The high luminosity and remarkable variability 
of CTD 135 and PKS 1413+135 observed by the Fermi/LAT, 
similar to blazars, demonstrate that their ?-rays are dominated by 
the radiation of a recently restarted aligned core-jet (Gan et al. 
2021, 2022). Theoretically, the mini-lobes of CSOs can also 
produce strong ^-ray emission via inverse-Compton (IC) 
processes (Stawarz et al. 2008; Migliori et al. 2014), or by 
thermal bremsstrahlung radiation (Kino et al. 2007, 2009), even 
through the hadronic processes (Kino & Asano 2011), similar to 
the large-scale lobes of Cen A and Fornax A (Abdo et al. 2010; 
McKinley et al. 2015). 

The core-dominance parameter, which is defined as the 
luminosity ratio (Rcg) of the core to extended regions, Rcg > 1 
for core-dominated and Rcg « 1 for lobe-dominated, is often 
used as an indicator of the beaming effect (e.g., Urry & 
Padovani 1995). We note that Cen A and Fornax A have low 
values of Rcg at the radio band compared with other y-ray 
emitting RGs. The Rcg values at 8 GHz’ are 0.48, 0.36, and 


^ The data are taken from the Radio Fundamental Catalog (RFC), http:// 


astrogeo.org/rfc/. 
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0.71 for PKS 1718-649, NGC 3894, and TXS 01284-554, 
respectively. Especially for PKS 1718-649, even no clear radio 
core is resolved (Angioni et al. 2019), resembling a mini- 
version of Fornax A (very weak radio core, McKinley et al. 
2015). Whereas, Rcg is ~1 for CTD 135 and 1.42 for PKS 
1413+135. It seems to be coincident with that PKS 1413+135 
and CTD 135 are reclassified as blazars rather than CSOs 
(Readhead et al. 2021; Frey et al. 2022). Different from CTD 
135 and PKS 1413+135, PKS 1718-649, NGC 3894, and TXS 
0128+554 definitely satisfy the CSO criteria asked by 
Kiehlmann et al. (2023). In this paper, we intend to explore 
the y-ray radiation properties of the three lobe-dominated CSOs 
by completely analyzing their Fermi/LAT observation data. 
Throughout, Hp — 71 km s! Mpc |, Qm = 0.27, and Q4 = 
0.73 are adopted in this paper. 


2. Lobe-dominated gamma-ray Emitting CSOs 


PKS 1718-649, also named NGC 6328 and located in a 
compact gas environment with redshift of z= 0.014 (Fosbury 
et al. 1977). It is also classified as a GHz-peak spectrum (GPS) 
source due to a turnover frequency of ~3 GHz at radio 
spectrum (Tingay et al. 1997). Its morphology at multiple radio 
bands is characterized as two sub-structures from southeast to 
northwest (Tingay et al. 1997, 2002), and the total projection 
size at 8.4GHz and 22.3 GHz is only 2.5 pc (Angioni et al. 
2019). No clear radio core is resolved at multiwavelength 
observations. However, the spectral index map from 8.4 to 
22.3 GHz (Figure 6 in Angioni et al. 2019) reveals that the 
region between two sub-structures displays a highly inverted 
spectrum and may be corresponding to the core, which cannot 
be observed at these radio frequencies due to synchrotron-self- 
absorption (Angioni et al. 2019). The estimated kinematic age 
is 70 + 30 yr via the apparent separation speed of 0.13 + 0.06c 
between two mini-lobes (Angioni et al. 2019, see also Giroletti 
& Polatidis 2009). 

NGC 3894, located in an elliptical galaxy with redshift of 
z — 0.01075 (Bilicki et al. 2014). Its radio spectrum peaks at a 
frequency near 5 GHz (Taylor et al. 1998). The high-resolution 
Very Long Baseline Array (VLBA) image at 1.4 GHz shows its 
radio core, the two-sided jets, and a southeast radio lobe with 
the total linear size of ~50 mas (Peck & Taylor 1998). The 
Very Long Baseline Interferometry (VLBI) observations at 
5 GHz covering 15 yr demonstrate that the two-sided pc-scale 
jets emerge from the radio core with a mildly relativistic speed 
of ~0.3c (Taylor et al. 1998). Principe et al. (2020) reported 
that its two-sided jets expand to southeast and northwest 
directions with an apparent separation velocity of 0.202 + 
0.005c between 1995 and 2017 and then estimated its 
kinematic age to be 59 + 5 yr. 

TXS 0128+554, hosted in an elliptical galaxy with redshift of 
z = 0.0365 (Huchra et al. 2012). The turnover frequency of its 
radio spectrum is less than 1 GHz (Lister et al. 2020). Its pc- 
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scale radio structure consists of an unresolved bright core and 
east—west side cambered lobes, and the overall projected extent 
of its lobes is similar in the multiple-radio-frequency images 
(Lister et al. 2020). With the observations at 2.3GHz and 
5 GHz, the faint steep-spectrum emission is revealed for the 
regions at 8.8 pc (eastern lobe) and 7.8 pc (western lobe) from 
the core (Lister et al. 2020). Lister et al. (2020) reported an 
apparent advance speed of 0.32+0.07c of its radio sub- 
structure and estimated its kinematic age of 82 + 17 yr under 
the assumption of a constant advance speed in the western lobe. 


3. Fermi/LAT Data Analysis and Results 
3.1. Data Analysis 


We use the public software fermitools> (ver. 2.0.8, Fermi 
Science Support Development Team 2019) with the binned 
likelihood analysis method to analyze the Fermi/LAT Pass 8 
data of the three CSOs. The data covering from 2008 August 4 
to 2022 May 1 (MJD 54682-59700) are extracted from the 
Fermi Science Support Center. We set the region of interest as 
a circle of radius 15? centered on the radio/optical positions of 
the three CSOs. The photon events in the energy range within 
0.1—300 GeV are considered. We filter the background y-ray 
contamination from the Earth limb by setting the maximum 
zenith angle as 90?. The good quality photon events are 
selected by a standard data quality selection criterion 
"(DATA QUAL»20)&&(LAT CONFIG—- 1)". The data are 
binned with a pixel size of 072 and 25 logarithmic energy bins. 
The instrument response function PSR3 SOURCE V3 is used 
in the data analysis. 

The Fermi/LAT 12yr Source Catalog (4FGL-DR3, 
Abdollahi et al. 2022), the diffuse Galactic interstellar emission 
(gll_iem_v07.fits), and the isotropic emission (iso_P8R3_SOUR- 
CE_V3_v1.txt) are included as background models. In the data 
analysis, we use the maximum test statistic (TS) to quantify the 
significance of the y-ray detection. TS is the logarithmic ratio of 
the likelihood of a model with the point source to that of a model 
without the point source, i.e., TS =2log(=*) (Mattox et al. 
1996). A threshold of TS > 25 is adopted to recognize a point 
source from the background. By subtracting all the background 
models, we first generate the 3° x 3° residual TS maps to 
investigate whether there are new y-ray sources not included in 
the 4FGL-DR3. No excess y-ray signal with TS > 25 is found in 
the residual TS maps, indicating that no new additional y-ray 
source is detected. We then produce the 2° x 2° TS maps of the 
three 4FGL point sources, 4FGL J1724.2-6501, 4FGL J1149.0 
+5924, and 4FGL J0131.2--5547. Using the gifindsrc tool, we 
also obtain the best-fit positions of the three 4FGL point sources, 
as shown in Figure |. Information on the three CSOs and the 
three associated 4FGL point sources is also given in Table 1. 


5. https: //fermi.gsfc.nasa.gov /ssc/data/analysis /software / 
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Figure 1. 2? x 2? TS maps of the three AFGL-DR3 point sources with the 68% and 9596 containments (black contours) of the best-fit positions. The cyan solid stars 
represent the radio/optical positions of the three CSOs. The maps are created with a pixel size of 0°05 and are smoothed. 


Table 1 
Fermi/LAT Analysis Results 
Source Name R.A.* Decl.* TS Lo.1.300 GeV D TS yar Association 
(deg) (deg) (erg s) 
PKS 1718-649 261.033 —65.0277 37.1 (1.02 + 0.19)E42 2.49 + 0.16 35.6 (~1.30) 4FGL J1724.2-6501 
NGC 3894 177.254 59.4158 96.4 (6.18 + 0.89)E41 2.16+0.11 37.3 (~1.50) 4FGL J1149.0+5924 
TXS 01284554 22.8277 55.7705 146.5 (1.62 + 0.17)E43 2.11 + 0.08 34.0 (1.10) 4FGL J0131.2--5547 


Note. 


è The best-fit position of the ~14 yr Fermi/LAT observation data obtained by the géfindsrc tool. 


During the spectral analysis, the normalization of the 
isotropic emission and the diffuse Galactic interstellar emis- 
sion, and the spectral parameters of the ^-ray sources listed in 
4FGL-DR3 within the range of 8° remain free, and other 
parameters are kept fixed. A spectrum should be significantly 
curved if TScury > 9 (corresponding to a 3c confidence level, 
Abdollahi et al. 2022), where TS curvy = 2 log( 2), Lip and 


Lp, are the hypothesis likelihoods of LogParabole model and 
power-law model testing, respectively. The TS.ury values of the 
three CSOs are much less than 9, therefore, the power-law 
spectral model is taken into account to fit their photon spectra 
in our analysis. The average spectra of the three CSOs with the 
power-law fits are presented in Figure 2, and the derived 
parameters are listed in Table 1. We also generate their long- 
term light curves with the power-law spectral fits, where the 
equal time interval of 180 days is adopted, as displayed in 
Figure 2. 

We also calculate the corresponding variability index TSya,, 
a common method to estimate the y-ray variability (Abdollahi 
et al. 2022). The definition of TS,,, is 


N 
TSyar = 2» /log(£;()) — logCCi(Fgoo))], (1) 
i=0 


where F; is the fitting flux for bin i, £;(F) is the likelihood 
corresponding to bin i, Foy is the best-fit flux for the global 
time by assuming a constant flux, and N is the number of time 


bins. We keep only the normalization of the target source free 
to vary and fix other parameters to the best-fit values during 
calculations. TS,4, = 53.0 corresponds to a 3c confidence level 
in the X JB distribution with N — 1 — 26 degrees of 
freedom, where N = 27 is the number of time bins. The y-ray 
emission of a source is probably variable when TS, > 53.0. 
The obtained TS,, values are listed in Table 1, and the 
corresponding confidence levels of variability are also given. 


3.2. Analysis Results 


PKS 1718-649. By analyzing the ~14yr Fermi/LAT 
observation data, a y-ray detection with TS ~37 is yielded 
for PKS 1718-649. As shown in Figure 1, the radio position of 
PKS 1718-649 falls within the 68% error circle of the best-fit 
position of 4FGL J1724.2—6501. Our result is coincident with 
the previous report that PKS 1718—649 is spatially associated 
with the y-ray point source 3FGL J1728.0—6446 (Migliori et al. 
2016), where 3FGL J1728.0—6446 is associated with 4FGL 
J1724.2-6501 in the 4FGL-DR3 (Abdollahi et al. 2022). PKS 
1718-649 has a steep spectrum at the GeV band with a photon 
spectral index of [',=2.49+0.16. Its ~l4yr average 
luminosity at the GeV band is (1.02 £0.19) x 10? erg s™!. 
No obvious variability is observed with TS,,, = 35.6, only four 
detection points are presented in the whole light curve, and the 
first three detection points are obtained with the first 7 yr 
observations. Using the 7 yr data set, Migliori et al. (2016) 
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Figure 2. Left column: the ~14 yr Fermi/LAT average spectra along the power-law function fit results. Right column: the long-term Fermi/LAT light curves with 
time bins of 180 days. The horizontal red dashed lines represent the ~14 yr average ^-ray luminosity of sources, and the red data points in the light curves denote the 
last residual time interval, only 157 days. The opened inverted triangles indicate the TS value of that energy/time bin less than nine. 


therefore found the y-ray signal of PKS 1718-649 with the 
same TS value (TS ~36) as ours. 

NGC 3894. A significant y-ray signal with TS = 96.4 is 
detected around NGC 3894. The radio position of NGC 3894 is 
located at the 68% containment of the best-fit position of 4FGL 
J1149.0+5924, as displayed in Figure 1. Its spectrum at the 
GeV band is flatter than that of PKS 1718-649, i.e., 
[,=2.16+0.11, which is slightly steeper than the derived 
value of I’, = 2.01 + 0.10 with the 10.8 yr observation data in 


Principe et al. (2020). NGC 3894 has the lowest y-ray 
luminosity among the three CSOs, being (6.18 + 0.89) x 
10?! erg s~". The confidence level of the flux variation is less 
than 2c, also indicating no variability in its y-ray emission. 
TXS 01284-554. The analysis with the ~14 yr Fermi/LAT 
observation data yields a significant y-ray detection with 
TS= 146.5 for TXS 0128+554. It has the highest ^-ray 
luminosity of (1.62 + 0.17) x 10? erg s^! and the flattest 
spectrum with I’, — 2.11 + 0.08 among the three CSOs. Same 
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Figure 3. The photon spectral index (IL) vs. the y-ray average luminosity (L,) 
in the LAT energy band. The data of blazars and RGs are taken from 4FGL- 
DR3. The data of CTD 135 and PKS 1413+135 are from Gan et al. 
(2021, 2022). 


as PKS 1718-649 and NGC 3894, no significant flux variation 
of y-rays more than the 2c confidence level is observed for 
TXS 01284-554. 

Note that Principe et al. (2021) investigated the ^-ray 
emission of some young radio sources by analyzing the 11 yr 
Fermi/LAT data. They reported that TS —36 and DL,— 
2.60 0.14 for PKS 1718-649, TS —95 and T,= 2.05 + 
0.09 for NGC 3894, TS = 178 and I’, = 2.20 + 0.07 for TXS 
0128--554, and found no evidence of significant ^-ray 
variability for the three CSOs. These results are consistent 
with our analysis results. 


3.3. The L, — I*, Plane 


As described above, the 4-ray average spectra of the three 
CSOs can be well fitted with a power-law function, none of 
them displays the significant flux variation in the GeV band, 
and all of them show the low ^-ray luminosity. These are 
dramatically different from the ^-ray emitting CSOs CTD 135 
and PKS 1413+135, both of them show the Log-parabola 
spectra, violent variability, and high luminosity in the GeV 
band, similar to blazars (Gan et al. 2021, 2022). The diversity 
of the y-ray emitting CSOs is also displayed in the L, — I, 
plane, where L, is the y-ray luminosity in the GeV band, and 
the data of blazars and RGs are taken from the 4FGL-DR3 
(Abdollahi et al. 2022). As illustrated in Figure 3, blazars 
together with CTD 135 and PKS 1413+135 occupy the high- 
luminosity area of the figure while RGs and three CSOs studied 
in this work are located at the low-luminosity region. It is 
widely known that blazars demonstrate frequent and violent 
variabilities, especially in the y-ray band, favoring a compact 
emission region with the significant Doppler boosting effect. 
Similarly, the high-luminosity y-ray emission of CTD 135 and 
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PKS 14134-135 is due to the episodic nuclear jet activities and 
is generated from a new-born aligned core-jet (Gan et al. 
2021, 2022). In contrast, RGs normally have the low 
luminosity and weak variability, and they are thought to be 
the parent populations of blazars with large viewing angles and 
small Doppler factors (Urry & Padovani 1995). Accordingly, it 
is still an open question whether the ^-rays of RGs originate 
from the jet with a large inclination angle or the large-scale jet 
extended regions. In the LAT band, Fornax A shows no 
significant variability while Cen A only displays some flux 
variation at a ^-3.5o confidence level in the low-energy band 
(Brown et al. 2017; Abdollahi et al. 2022). In particular, the 
Fermi/LAT observations definitively demonstrate the ^-rays 
from the large-scale lobes of Cen A and Fornax A, to account 
for more than 50% and 86% of their total y-ray flux, 
respectively (Abdo et al. 2010; Ackermann et al. 2016). 
Thereby, the 4-rays of Cen A and Fornax A are lobe- 
dominated. Recently, the ^-ray emission from FR II RG IGR 
J18249-3243 was also found to be lobe-dominated (Bruni et al. 
2022). In view of no obvious ^-ray variability and the lobe- 
dominated radio emission for TXS 0128+554, PKS 1718-649, 
and NGC 3894, similar to these RGs, we speculate that their y- 
rays stem from the extended mini-lobes. 


4. Discussion 
4.1. Broadband SED Modeling and ^j-ray Origin 


In order to further explore the y-ray emission properties of 
the three CSOs, we collect the low-energy band data from 
literature and the NASA/IPAC Extragalactic Database (NED) 
as well as the ASI Science Data Center (ASDC) and construct 
their broadband spectral energy distributions (SEDs), as 
displayed in Figure 4. A detailed description of data is in the 
caption of Figure 4. Considering the similar radiation properties 
of the three CSOs to RGs Cen A and Fornax A, the lobe- 
dominated emission in the radio band, and no significant 
variability in the GeV band, we thus infer that the 4-rays of the 
three CSOs are also dominated by the extended region. Due to 
the broad second bump in the broadband SEDs of the three 
CSOs, we thus intend to represent their SEDs with the two- 
zone leptonic model, similar to that in other y-ray emitting 
compact radio sources (e.g., Zhang et al. 2020; Gan et al. 2021, 
2022; Gu et al. 2022). The synchrotron (syn), synchrotron-self- 
Compton (SSC), and external Compton (EC) scattering of the 
relativistic electrons in both the core and extended region are 
considered. The synchrotron-self-absorption (SSA) effect at the 
low-energy band, the Klein-Nishina effect at the high-energy 
band, and the absorption of high-energy y-ray photons by 
extragalactic background light (EBL; Franceschini et al. 2008) 
are also taken into account during the SED modeling. 

The electron distributions in both the core and extended 
regions are assumed as a broken power-law, which are 
characterized by a density parameter No, a break energy %, 
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Figure 4. Observed SEDs with the two-zone (“1” for the extended region and “2” for the core region) leptonic model fitting results. The data marked as opened gray 
squares are taken from the NED, ASDC, and literature (Balasubramaniam et al. 2021; Sobolewska et al. 2022). The red squares and inverted triangles denote the 
Fermi/LAT average spectra obtained in this work. The olive stars (core) and the opened cyan circles (lobe) represent the flux of core and lobe, respectively, which are 
collected from the NED and literature (Tingay et al. 1997; Taylor et al. 1998; Lister et al. 2020). The opened blue squares indicate total flux taken from Lister et al. 
(2020). The blue dashed-dotted lines show the extrapolated VHE spectra from the Fermi/LAT spectra, while the orange dashed-dotted line in TXS 0128+554 panel 
represents the extrapolating VHE spectrum considered the EBL absorption. The data of the soft X-ray band between two vertical lines for NGC 3894are not 


considered during SED modeling. 


and indices p, and p; in the energy range of [^i Ymax] The 
radiation region is assumed as a sphere with radius R, magnetic 
field B, and Doppler factor 6. Since the mildly relativistic 
motions have been observed and reported for NGC 3894 
(Principe et al. 2020) and TXS 01284-554 (Lister et al. 2020), 
we take their relativistic effect into account for both the core 
and extended regions, i.e., ó — 1.3 and I = 1.04 for NGC 3894 
(Principe et al. 2020), 6 = 1.2 and I = 1.06 for TXS 0128+554 
(corresponding to a viewing angle of 0 = 52°, Lister et al. 
2020), where I is the bulk Lorenz factor. The jet is accelerated 
and collimated reaching up to sub-pc/pc scale (see Boccardi 
et al. 2017 for a review). After the acceleration phase, the jet 
assumes a conical shape and initiates deceleration as a 
logarithmic function of distance (r) from the black hole (see 
Potter & Cotter 2013; Wang et al. 2022). We assume the same 
value of I' for both the core and extended regions due to the 
total projection size of NGC 3894 and TXS 0128+554 being 
less than 20 pc. Assuming the same advance speed for the two- 
sided mini-lobes, an apparent advance velocity of only ~0.06c 
is observed for PKS 1718-649 (Angioni et al. 2019), we thus 
take I’ = ô= 1 for both the core and extended regions during 
SED modeling. For the core region, we take R = 6cAt/(1 + z) 
and At= 180 days. Due to the absence of apparent variability 
in the ^-ray light curves, the time-bin value of the light curves 
is adopted, similar to other compact radio sources emitting ~- 
rays (Zhang et al. 2020). In this case, R » 0.1 pc, the energy 
dissipation region should be located at ~10 * R from the black 
hole and outside of the broad-line regions, and the photon field 
of torus provides the dominant seed photons of the EC process 
(see Ghisellini & Tavecchio 2009). The overall projection size 
of the radio morphology of the three CSOs is all less 
than 20 pc. The seed photon field of the EC process is still 
dominated by the photons from torus within this scale 


(Ghisellini & Tavecchio 2009, see also Stawarz et al. 2008). 
We thus only consider the external photon field of torus to 
calculate the EC process for both the core and extended region. 
The typical energy density of the torus photon field 
of Ug —3 x 10 * erg cm ? (Cleary et al. 2007; Kang et al. 
2014) is taken for the core region, whereas the values of Ur for 
the extended regions are estimated according to Equation (21) 
in Stawarz et al. (2008) and are given in Table 4. For the core 
region, 4,,, = 1 and 4,,, = 10^ are set. p; is roughly estimated 
by the spectral index of the X-ray observation data while p; — 4 
is fixed. W, No, and B are free parameters. For the extended 
region, R is equal to the total projection size of the radio 
morphology, as listed in Table 2. ^;.. is taken as 300 or 500 to 
match the observation data at the radio band and ^,,,, is set as 
Ymax = 100^, (roughly constrained by the last upper-limit data 
point in the Fermi/LAT spectrum). p, and p» are constrained 
by the radio and y-ray spectral indices, respectively. y», No, and 
B are also free parameters. As shown in Figure 4, the 
observational data at the near-IR-optical band are dominated 
by the starlight of host galaxy. The data of the soft X-ray band 
(0.5—1.5 keV) for NGC 3894 are not corrected the absorption 
(Balasubramaniam et al. 2021). We thus do not take these data 
into account during the SED modeling. Due to the limited 
observation data, we cannot give a constraint on the fitting 
parameters, and we only obtain a set of model parameters that 
can make an acceptable fit. We also plot a cartoon illustration 
for the radiation model, as depicted in Figure 5. 

The fitting results are presented in Figure 4 and the modeling 
parameters are given in Table 2. The broadband SEDs of the 
three CSOs can be well explained by the two-zone leptonic 
radiation model. Under this scenario, the y-rays are produced 
by the IC process of the relativistic electrons in the extended 
region and the X-rays are dominated by the radiation of the 


Table 2 
SED Fitting Parameters 
Compact Core Extended Region 
Source R B 6 r 2m Yb Yoi No Di pa R B ô DL frais Yb ax No Di P2 
(cm) (mO) (em ?) (Pc) (mG) (em 7) 
PKS 1718-649 4.60E17 80 1 1 1 1706 — 1.0E4 1321 2.0 4.00 2.5 42 1 1 300  9.5E3  9.5E5 50 2.33 4.00 
NGC 3894 6.00E17 10 13 1.04 1 1450 1.0E4 6.26E-2 1.08 4.00 10.9 1.7 13 1.04 500 5.0E3 5.0E5 5.5E-6 1.02 3.20 


TXS 0128+554 — 5.40E17 10 12 1.06 1 3056  1.0E4 477 2.00 4.00 17.9 1.0 12 1.06 500 2.5E4 25E6 7.5E-5 1.25 3.25 
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Figure 5. A cartoon illustration of the two-zone leptonic radiation model. The distance between the core region and the black hole is conservatively estimated to be 
greater than 1 pc, with a rough approximation of 10 * R. The distance between the extended region and the black hole is approximately half of the total projection size 
of the radio morphology, where the term “LS” refers to the overall projection size of sources. 


core region. As listed in Table 2, the derived B values of core 
regions are higher than that of extended regions while ^ of the 
core is smaller than that of the extended region. We find that 
the magnetic field strengths between the core region (B4) and 
the extended region (B2) for the three CSOs basically satisfy 
the relation of B4/B; = rz/r,, where r, and r» are the distances 
from the core and the extended regions to the black hole, 
respectively. This result is consistent with the VLBA observa- 
tions for some AGN jets (O'Sullivan & Gabuzda 2009). The 
derived B values of extended regions are several mG, being in 
agreement with the typical values measured in CSOs (Orienti 
et al. 2006; Lister et al. 2020). As listed in Table 2, ^i... ~ 106 
is required to account for the y-ray emission originating from 
these extended regions of CSOs. The detection of GeV ^-rays 
in the large-scale lobes of RGs Cen A and Fornax A (Abdo 
et al. 2010; Ackermann et al. 2016) provides compelling 
evidence for the acceleration of ultrarelativistic electrons in the 
extended regions of AGN jets, particularly highlighted by the 
very high-energy detection along the jet of Cen A (H.E.S.S. 
Collaboration et al. 2020). The leptonic radiation model 
employed to explain that the observed y-ray flux from these 
extended regions of AGN jets requires electron Lorentz factors 
on the order of 10° (Persic & Rephaeli 2019), or even higher 
(107-105; H.E.S.S. Collaboration et al. 2020). 

It should be noted that we used a phenomenological electron 
spectrum, a broken power-law spectrum, to model the SEDs 
here. Although a power-law electron spectrum with an 
exponential cut-off would be more naturally motivated by the 
particle acceleration theory, a broken power-law electron 


spectrum can better represent the symmetric peaks in the 
observed SEDs (e.g., Aharonian et al. 2009; Ghisellini & 
Tavecchio 2009; Chen et al. 2012; Zhang et al. 2012; Aleksié 
et al. 2014). It is unclear whether the electron distribution has a 
low-energy cutoff with Yain > 1 in theory; generally, it is fixed 
as Yuin = 1 (e.g. Ghisellini & Tavecchio 2009; Zhang et al. 
2012). Sometimes, 7,,;, can be roughly constrained with X-ray 
data (e.g., Tavecchio et al. 2000; Zhang et al. 2014), or a low- 
energy cutoff in the electron spectrum is required to explain the 
radio spectra of some large-scale jet substructures, with ^... 
being several hundreds (Harris et al. 2000; Hardcastle et al. 
2001; Lazio et al. 2006; Godfrey et al. 2009). During SED 
modeling, 7y,,,, normally is constrained by the last observation 
point in the 4-ray band. The theoretical upper limit of the 
maximum electron energy 7,,,, can be determined by tacc = 
min (fac, f;59)), Where tacc, teses tooo! are the timescales of 
acceleration, escape and cooling, respectively. We estimated 
the three theoretical timescales using the Equations from Wang 
et al. (2022), and found that the derived values of ^... from 
SED modeling are much smaller than the theoretical upper 
limit of Yax: In addition, % derived from SED modeling is not 
exactly the same as the cooling break since ^ results from a 
complex physical process, including the adiabatic losses, the 
particle escape, and the cooling (e.g., Ghisellini & Tavecchio 
2009). 

The energy densities of the non-thermal electrons (Ue) and 
magnetic fields (Ug) of extended regions are also listed in 
Table 4. The ratio of U./Ug of the three CSOs is close to that 
of the lobes in Cen A and Fornax A, which are also derived by 
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Table 3 
Derived Parameters of the Core Region 
Source Us U: U:/Ug Ps P. P, Pe 
(erg cm ?) (erg cm ?) (erg s7!) (erg s^) (erg s (erg s ! 
PKS 1718-649 2.55E-4 2.59E-3 10.2 5.07E42 5.17E43 8.72E41 5.76E43 
NGC 3894 3.98E-6 6.54E-5 16.4 1.46E41 2.40E42 6.15E41 3.16E42 
TXS 01284554 3.98E-6 3.31E-3 832.1 1.23E41 1.02E44 1.18E43 1.14E44 
Table 4 
Derived Parameters of the Extended Region 
= 
Source VIR Ur Ug Ue Ue/Ug Pg Pe P, Pret 
(Hz) (erg cm ?) (erg cm ?) (erg cm ?) (erg s!) (erg sl) (erg s7!) (erg s t 
PKS 1718-649 2.54E13 6.99E-7 7.02E-7 1.38E-5 19.7 3.93E42 7.76E43 2.11E42 8.36E43 
NGC 3894 2.97E13 5.83E-8 1.15E-7 3.31E-8 0.3 1.33E43 3.82E42 1.44E41 1.72E43 
TXS 0128+554 2.59E13 9.39E-9 3.98E-8 2.51E-7 6.3 1.28E43 8.12E43 6.07E42 1.00E44 


the SED modeling (Persic & Rephaeli 2019), and is also in 
agreement with the values in the radio lobes of some RGs, 
which are obtained by the X-ray measurements (e.g., Croston 
et al. 2005; Kataoka & Stawarz 2005; Isobe et al. 2011). 
Based on the SED fitting parameters, we estimate the jet 
powers (Pe) of the core and extended regions in the case of the 
electronpositron pair jets. The powers of the non-thermal 
electrons (Pe) and magnetic fields (Pg) are calculated with 
P= nR'T?cU;, where U; can be U. or Ug. The radiation power 
of jets is estimated by P, = R?I?cU, = Lust / 46^, where Loot 
is the total luminosity of non-thermal radiation from a radiation 
region. pr. is the sum of Pe, Pg, and P,. As given in Tables 3 


and 4, Pe of the core is within 3 x 10*°-104 erg s^! while P 
of extended regions are in the range of ~10*-10“ erg s !.Itis 
found that the jet powers of the core and extended regions for 
the three CSOs are roughly of the same order, and similar 
results have been reported for a y-ray-emitting RG 3C 120 
(Zargaryan et al. 2017) and a blazar 4C +49.22 (Zhang et al. 
2018). We also estimate the kinetic power (Pķkin) of the three 
CSOs with the radio luminosity at 1.4 GHz of extended regions 
using Equation (1) in Cavagnolo et al. (2010) and obtain 
1.4 x 10“ erg s^!, 3.0 x 10? erg s^! and 7.3 x 10? erg s^! 
for PKS 1718-649, NGC 3894 and TXS 01284-554, 
respectively. Py;, is also roughly of the same order as Pe, of 
the core and extended regions. 


4.2. Possible VHE Emission 


Recently, CSO PKS 1413+135 is detected at the very high 
energy (VHE) band by the MAGIC telescope (Blanch et al. 2022; 
Gan et al. 2022). Although the y-rays of PKS 1413+135 should 
be dominated by an aligned jet emission and may be totally 
different from the three lobe-dominated CSOs, we are still 
wondering whether the three CSOs could be detected at the VHE 


band. In particular, the H.E.S.S. Collaboration reported that Cen 
A has the VHE emission along its large-scale jet, providing direct 
evidence for the presence of high-energy electrons in the extended 
regions of AGN jets (H.E.S.S. Collaboration et al. 2020). 

Extrapolating the Fermi/LAT ^-ray spectrum in the VHE 
band is a common method to estimate the intrinsic VHE 
spectrum of sources (Paiano et al. 2021; Zhu et al. 2021; Malik 
et al. 2022). As described in Section 3.2, the spectra observed 
by the Fermi/LAT for the three CSOs are well represented by a 
power-law function, and no evidence of variability has been 
found in their long-term ^-ray light curves. We therefore 
directly extend the power-law spectra at the GeV band to the 
VHE band for the three CSOs, as illustrated in Figure 4. The 
sensitivity curves of Cherenkov Telescope Array south and 
north (CTA-S and CTA-N, 50 hr) are also presented in 
Figure 4, where the sensitivity curves are obtained from the 
CTA webpage.^ It can be found that the extrapolated flux at the 
VHE band of PKS 1718-649 and NGC 3894 is below the 
sensitivity curve of CTA. Especially for PKS 1718-649, 
together with its steep spectrum at the GeV band, it would be 
difficult to be detected by the CTA in the future. The 
extrapolated intrinsic VHE spectrum of TXS 0128--554 clearly 
surpasses the sensitivity curve of the CTA; however, when 
considering the EBL absorption, it only marginally intersects 
with the CTA’s sensitivity curve, as depicted in Figure 4. 
Nevertheless, theoretical suggestions have proposed episodic 
jet activity for compact radio sources (e.g., Reynolds & 
Begelman 1997; Czerny et al. 2009), and the episodic nuclear 
jet activity has been reported for TXS 0128+554 (Lister et al. 
2020). The reactivation of TXS 0128+554 could position it as 
a potential candidate for VHE emission, making it detectable 
by the CTA. 


6 https: //www.cta-observatory.org /science /cta-performance / 
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5. Summary 


We completely analyzed the ~14 yr Fermi/LAT observation 
data of three radio-lobe-dominated CSOs (PKS 1718-649, NGC 
3894, TXS 0128+554) to explore their y-ray emission 
properties. They are all characterized by a low y-ray luminosity 
(10*'-10* erg s^, their average spectra in 0.1—300 GeV can be 
well fitted by a power-law function, and no significant variability 
is observed in their long-term ^-ray light curves. These ^-ray 
features are prominently different from those of the ^-ray 
emitting CSOs PKS 1413+135 and CTD 135. In the L, — I, 
plane, the three CSOs are located at the region occupied by RGs, 
also distinctly different from PKS 14134-135 and CTD 135, 
which are similar to blazars. Considering the similar radiation 
properties between the three CSOs and the ^-ray emitting RGs, 
specifically similar to RGs Cen A and Fornax A, together with no 
significant variability in the y-ray band and no superluminal 
motion at the radio band, we speculated that their ^-rays stem 
from the lobe-dominated radiation, different from CTD 135 and 
PKS 14134-135 (Gan et al. 2021, 2022). Their broadband SEDs 
can be represented by the two-zone leptonic model and 4-rays are 
produced by the IC processes of the relativistic electrons in 
extended regions. The ^-rays are produced by the IC process of 
the relativistic electrons in extended regions and the X-rays are 
dominated by the radiation of the core region. This radiation 
model can be checked by studying whether the flux variations in 
X-ray and ^-ray bands are correlated. Based on the derived fitting 
parameters, we calculated the jet powers of the core and extended 
regions by assuming the electron-positron pair jet and found that 
the jet powers of the core and extended regions for the three 
CSOs are roughly of the same order. TXS 0128+554 may be a 
VHE emission candidate and be detected by CTA in the future. 
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